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The particle-hole excitation spectrum for doped graphene is calculated from the dynamical po- 
larizability. We study the zero and finite magnetic field cases and compare them to the standard 
two-dimensional electron gas. The effects of electron-electron interaction are included within the 
random phase approximation. From the obtained polarizability, we study the screening effects 
and the collective excitations (plasmon, magneto-excitons, upper- hybrid mode and linear magneto- 
plasmons). We stress the differences with the usual 2DEG. 



I. INTRODUCTION 

The particle-hole excitation spectrum (PHES) of a ma- 
terial is an extremely useful ingredient for the under- 
standing of its electronic properties, namely at low ener- 
gies. It reveals, for example, collective modes such as the 
plasmon, the electric polarizability of the material, and 
the screening properties of the electrons in it. In this pa- 
per, we review the particular PHES of two-dimensional 
electron systems, both for non-relativistic electrons in 
a usual two-dimensional electron gas (2DEG)^ and for 
massless electrons in graphene. 2 The main focus of the 
present review article is on the PHES and the collective 
modes of doped graphene in a strong perpendicular mag- 
netic field in the integer quantum Hall regimeP 4 * 5 * 6 * 7 * 8 * 9 * 
For pedagogical reasons, we compare them to the cor- 
responding results of the 2DEG 1011 as well as to the 
PHES in graphene without a magnetic fiel d, which h as 
been theoretically studied in detail recently!- 12 1 13 1 14 1 15 1 

The main difference between the PHES in graphene 
and that in the 2DEG, in a strong magnetic field stems 
from the quantization of the electrons' kinetic energy in 
Landau levels (LL): these levels are equidistantly spaced 
in the case of non-relativistic electrons in the 2DEG, 
e 2DEG = oj c (n + 1/2), in terms of the cyclotron frequency 
co c = eB/mb, where is the band mass (we use a sys- 
tem of units such that h = 1). In contrast to this, the LLs 
in graphene occur in two copies as a consequence of the 
two energy bands (A = + for the conduction and A = — 
for the valence band), and their level spacing decreases 
with increasing LL quantum number n, 

Vp i — / 

e\^ n = Xe n = X-—v2n oc WBn . (1) 

where is the Fermi velocity in graphene and Ib = 
1/y/eB is the magnetic length. This difference in LL 
quantization as well as the absence of backscattering due 
to the chirality properties of electrons in graphene^ yield 
a strikingly different PHES for graphene when compared 
to the 2DEG. Whereas, in the latter, the collective exci- 
tations are dominated by essentially horizontal weakly- 
dispersing magneto-excitons pl in addition to the up- 



per hybrid mode, those in graphene are linear magneto- 
plasmons 8 that disperse roughly parallel to an energy line 
uj = Vj?q, as a function of the wave vector q. The pre- 
cursors of these modes are already visible in the PHES 
for non-interacting electrons and acquire coherence once 
electron-electron interactions are taken into account, e.g. 
on the level of the random-phase approximation (RPA). 

The paper is organized as follows. In Sec. II, we intro- 
duce the basic expressions for the polarization function 
of graphene in a strong magnetic field. The intermediate 
steps of the derivation may be found in Appendix A. Sec- 
tion III is devoted to a discussion of the PHES in the stan- 
dard 2DEG for non-relativistic electrons, whereas that 
for doped graphene is presented in Sec. IV. Both sec- 
tions comprise also a short review of the B = case. 
In Sec. V, we aim at a physical interpretation of main 
features of the two different PHES in a strong magnetic 
field within a wave-function analysis, and the screening 
properties in the static limit are reviewed in Sec. VI. 

II. POLARIZABILITY 

The Hamiltonian for graphene in a magnetic field can 
be expressed as Hq — Vp7r • a for the K valley and 
Hq' = v f tt - a* for the K' valley, where a = (a x , cr y ) 
are Pauli matrices, and 7r = p + eA is the gauge- 
invariant momentum with p = — iV, and A is the vec- 
tor potential. In the symmetric gauge, the latter reads 
A = (— By/2 : Bx/2, 0), where B is the modulus of the 
magnetic field that we choose in the ^-direction. The 
Fermi velocity v F = 3ta/2 is expressed in terms of the 
nearest-neighbor hopping integral t ~ 3 eV and the 
carbon-carbon distance A. The eigenstates com- 

ponents |^Ca;A,n,£) are: 

|^+A;A,n,£) = \^-B-\,n,l) = ~^l*|n - l,t) 
\^-A-X,n/) = \^B-X,n,l)=2* n \nJ) (2) 

where n is a positive integer, A = ± for states of pos- 
itive/negative energy and A = for n = 0. The cor- 
responding eigenenergies are given in Eq. 0. The in- 
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dex ( = +(— ) denotes electrons in the K(K / ) valleys 
and a = A(B) the A(B) sublattice component of the 
electronic wave function. We have furthermore intro- 
duced the simplified notation 1* = y/(l — 5 n> o)/2 and 
2* = + 5 n $)/2. Here the quantum number n labels 
the LL, whereas the other quantum number which de- 
termines the LL degeneracy, varies from to Nb — 1, 
with Nb = Atib = A/2nl 2 Bl in terms of the total sample 
surface A . The states \n,£) are the eigenvectors of the 
Hamiltonian 



a/2DEG 

n 



2m6 



(3) 



for the standard 2DEG in a magnetic field. 
The bare polarization function 



iII (q, W ) 



duj'dk 
"(2^)3" 



Tr [G°(k,w')G°(k + q,w'- 



may be calculated with the help of the Green's func tions 
G°(k, uS) for non-interacting electrons [see Eq. (A5)] and 
reads, for the case of a strong magnetic fieldP 



n°(q,-) = EE 



e(A'c„0 - e(A£ n ) ^av 

:^ nn'\H)' \°) 



XX' n,n 



< t \i n - A'f n / + U + iS 



-XX' , 



The expression of the functions T nn , (q) and the details of 
the calculation may be found in AppendixjA] II contains 
two separate contributions, 



N F 



n° (q, W ) = £ (q, u>) + IF- (q, u>) . 



(6) 



The vacuum contribution LP ac (q, uo) takes into account 

inter-band processes, whereas Y^n=i ^n F (Q? °°) repre- 
sents intra-band contributions when the Fermi energy 
€f — £n f lies in the conduction band, as we assume 
implicitly from now on. 



A. Effect of electron-electron interaction 

From II (q, uo) we can calculate the renormalized po- 
larization function in the RPA, which is defined as 

u RPA (a , _ n°(q,a>) II°(q >tt ,) 

£ RPA {%L0 ) l- V (q)n°(q, W ) {> 

where v(q) is the unscreened two-dimensional Coulomb 
potential in momentum space 
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^|q| 



(8) 



in terms of the background dielectric constant and 
e RPA (c±,uj) = 1 — u(q)II (q, uj) is the dielectric function. 



Long-range electron-electron interaction usually leads to 
the appearance of collective modes in the spectrum, such 
as the plasmon, the dispersion of which is defined from 
the zeros of the dielectric function 



(q,") = 0. 



(9) 



Collective modes will be discussed in detail for both, 
graphene and a standard 2DEG, in the following sections. 



III. PARTICLE-HOLE EXCITATION 
SPECTRUM OF A STANDARD 2DEG 

In this section we briefly review the results for the 
PHES in a 2DEG and start with the case of zero magnetic 
field. The polarization function Eq. Q for a system of 
free electrons with parabolic band £k = k 2 /2mi ) can be 
expressed, after some manipulation, 

aPl 



n 2DEG(q^) = p( e F) 



■iS 



v F q 



v F q 



Q_ 

2k F 

2k F 



(10) 



where p(ep) is the density of states at the Fermi level, 
S ~ r _1 accounts for a finite life time o f the quasiparti- 
cles, and *&(z) = z — sgn(Re z)\J z 2 — 1. In Fig. [lja) we 
show a density plot of ImIl2 DEG . From the imaginary 
part of n° we can obtain the PHES, which is the region 
of the momentum-energy plane where it is possible to 
excite electron-hole pairs. This region is defined, for a 
non-interacting electron gas in the a bsence of a magnetic 
field, as the continuum of Fig. 1(a) with non-zero Imll , 
which corresponds (for S — > 0) to the region delimited 
by the solid black lines. The boundaries of the spec- 
trum are defined by max[0, uo- {q)\ < \uo\ < u;+(q), where 
uj±(q) — q 2 /2ml) ± v F q, in terms of the Fermi velocity 
v F = de/dk\k F = kp/mb — \/2^Jm\ ) ^ which in contrast 
to graphene depends on cf- Notice that this PHES is 
not uniform, but presents some structure. Apart from 
the quasi-homogeneous yellow region, two other zones 
are worth describing: the blue one, with a strong spec- 
tral weight, which is the precursor of the plasmon mode, 
as we will see below, and the almost red low energy re- 
gion, with a very weak spectral weight and which, in a 
ID system, would belong to the forbidden zone of the 
spectrum for particle-hole excitations. 

Electron-electron interactions lead to the appearance 
of a collective mode, the plasmon, which can be captured 
within the RPA. The imaginary part of U RPA (q, a;), 
shown in Fig. 1(b), reveals above the continuum bound- 
aries a well-defined peak centered at uo = uo pi the fre- 
quency of the plasmon mode. The broadening of the 
peak depends on the value of 5 ~ r _1 and is due, e. g., 
to scattering of electrons by disorder. The dispersion re- 
lation of the plasmon may be calculated from the zeros 
of the RPA polarization function. From the q — > and 
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FIG. 1: Density plot of ImII(q, uS). Plots (a) corresponds to 
non-interacting polarization of a 2DEG, Eq. (10), whereas 
(b) accounts also for electron-electron interactions in the 
RPA. Plots (c) and (d) show the imaginary part of the non- 
interacting and RPA polarization function, respectively, of 
a 2DEG in a magnetic field. In (a) and (b) we have used 
6 = 0.2e F . In (c) and (d), Nf = 3 and S = 0.2oo c . In (b) and 
(d), we have used r s c± 3. 



uo <C v^q expansion of the polarization function, the long 
wavelength limit of the plasmon dispersion is found to 



w p (q) 



l2e 2 e F 



(11) 



where the Fermi energy in a 2DEG with a parabolic band 
can be expressed in terms of the the uniform density of 
electrons n e \ as ep = 7rn e //m6. Furthermore, from the 
analytic solution of Eq. ([9| , an exact dispersion relation 
of the plasmon mode may be obtained to all orders in 
q (see Ref. [T8|) . At low energies, it disperses as y/q, 
and at some critical wave vector q c , the mode touches 
the electron-hole continuum. This critical value may be 
obtained from 18 



(gc/fcf) 2 , {q c /k F f 
V2r s 



4r? 



= 1, 



(12) 



in terms of the dimensionless interaction parameter r s = 
2m\ ) e 2 1 'e^kp. Above q c , the plasmon disperses roughly 
parallel to the boundary of the particle-hole contin- 
uum while being Landau-damped due to its decay into 



electron-hole pairs [see Fig. 1(b) 



In the presence of a strong magnetic field perpendicular 
to the 2DEG, the bare polarizability can be expressed 



ttO / \ V^ 7 J~ n,m(q) , / + -\ 



m=l 
•\Nf — 1 



uj — rnoj c + iS 
and uo + 



(13) 

where Yl = En=max(o,Ar F -m) and ^ + ^ indicates 
the replacement uo + iS — > —uo — iS. The form factor due 
to the wave-function overlap is 



3~ n,m(Q) — 



2tt/| 



(n - 



m) 



2/2 
B 



qH 



a 2 ] 2 
q t B 



(14) 

A density plot of Imll is shown in Fig. 1(c) for Np = 3. 
In the presence of a strong magnetic field, ImII (q, uj) is 
a sum of Lorentzian peaks centered at uo = mw c , with 
m = n 1 — n > 1, the difference between the LL indices 
of the electron n 1 and the hole n. Therefore the PHES is 
chopped into horizontal lines, separated by a constant en- 
ergy uo c . The width of each horizontal line is proportional 
to the disorder broadening of the Landau levels. (The 
peaks become (5-functions in the clean limit 5 — » + .) No- 
tice that within each horizontal line, the spectral weight 
is not homogeneously distributed; one observes inde ed a 
superstructure of N F + 1 brighter regions in Fig. |l(c)| 
following lines parallel to the edges of what used to be 
the particle-hole continuum in the zero-field case. These 
edges delimit, also for non-zero values of the magnetic 
field, the PHES region of non- vanishing spectral weight. 

In the presence of electron-electron interactions, the 
density fluctuation spectrum of a 2DEG at integer filling 
factors is dominated by a single set of collective modes, 
known as horizontal magnet o-excitons. The frequency of 
these modes tend to muo c (where m is a positive inte- 
ger) in the long wavelength limit, and their dispersion 
was calculated by Kallin and Halperin in the RPA and 
in the time-dependent Hartree-Fock approximation El In 
Fig. 1(d) we show the RPA excitation spectrum for 
r s = 2m\ ) e /sbkp ~ 3. Furthermore, one notices that the 
plasmon energy is renormalized by the magnetic field and 
evolves into the so called upper hybrid mode. Its disper- 
sion relation can be expressed as^ 



(15) 



where an approximate expression for uo p (q) has been 
given in Eq. (11 ). 



IV. PARTICLE-HOLE EXCITATION 
SPECTRUM OF DOPED GRAPHENE 

We start with a review of the polarization function in 
the absence of a quantizing magnetic field, which turns 
out to be useful in the understanding of the PHES also 
at B 7^ 0. In the continuum approximation, the non- 
interacting polarization function at zero temperature can 
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FIG. 2: Same as Fig. ^ but for the case of graphene. The 
solid lines delimitate the boundaries of the PHES. Region 
I corresponds to the intra-band and region II to the inter- 
band region of the PHES of doped graphene. In (a) and (b), 
S = O.Isf- In (c) and (d), 5 = 0.2vf/Ib- In (b) and (d), we 
have used r s — 1. 



be calculated from 

g e(A'& +q ) - e(A&) 



n°( q ,w) 



iS 



F AA ,(k,k+q), 
(16) 

where A£k = Avp|k| — ep is the quasiparticle energy rel- 
ative to the Fermi level, g — g s g v = 4 accounts for spin 
and valley degeneracy and 



F AA ,(k,k + q) 



1 + AA' COS # k ,k+q 



(17) 



is the chirality factor or wave function overlap, where 
#k,k+q is the angle between k and k + q. n°(q, uo) 
was calculated for undoped graphene in Ref. [2Ql where 
it was found that ImII (q,cc;) ~ q 2 (uo 2 - Vp<? 2 )~ 1/2 , 
which implies that the massless non-interacting electrons 
in graphene have an infinite response at the threshold 
uo = v^q. The reason for this behavior is twofold: first 
the threshold is determined by the linear dispersion re- 
lation and second, the chirality factor (17) suppresses 



backscattering in graphene. This feature is still present 
in doped graphene, although the form of the pola rization 
function is richer than in the absence of doping,^ 2 * 14 * 15 * 
as may be seen in Fig. [2j where we show a density plot 
of II(q, uo) for doped graphene at B = 0. One notices 
that most of the spectral weight is actually concentrated 
around uo = v^q, as one would expect from the suppres- 
sion of backscattering (2k p processes). This is similar to 
the case of zero doping. 



In doped graphene, however, there are two regions 
of non-vanishing spectral weight which arise from intra- 
band (region I) and inter-band processes (region II). 
These two regions are separated by the diagonal line uo = 
v^q. For zero doping, there are naturally only inter-band 
processes. The intra-band contributions are restricted to 
region I the boundaries of which are max(0, v^q — 2ep) < 
\uo\ < v^q. This is the only kind of processes present in 
the 2DEG, discussed in the preceding section. The fea- 
tures of this zone are, apart from the different shape of 
the boundaries in the two cases, similar to the PHES of 
the 2DEG, although the spectral weight is no longer ho- 
mogeneously distributed within this region but vanishes 
when approaching the right boundary as a consequence of 
the chirality factor. Inter-band particle- hole excitations 
are restricted to region II in the PHES, with boundaries 
| a; | > max(v F 3S —vfQ + 2ep). Indeed, the inter-band ex- 
citation of lowest energy needs to overcome the Fermi 
energy, which is therefore the lower bound of region II 
and which approaches zero in undoped graphene, where 
the region of inter-band excitations covers the whole part 
uo > vpq of the spectrum. Furthermore, due to the pres- 
ence of two energy bands, there is the possibility of direct 
transitions from the valence to the conduction band with 
q = momentum transfer. The q = transition of low- 
est energy involves an energy cost that is twice the Fermi 
energy, which yields a gapped region in the PHES in Fig. 
|2|a), denned by < uo < 2s F . Notice however that the 
q = transition is suppressed due to the chirality factor 



(17). 



The real and imaginary parts of the polarization func- 
tion are related via the Kramers-Kronig relations. Their 
behavior is shown in Fig. [3] where we plot — H°(qo,uo) 
for different wave- vectors qo for both, 2DEG [plots (a) 
and (b)] and graphene [plots (c) and (d)]. We note that 
most of the spectral weight (region with non-zero Imll ) 
of graphene is concentrated near uo = v-pq (due to the ab- 
sence of backscattering) , whereas in the 2DEG it is more 
uniformly distributed. 

When electron-electron interactions are included in the 
problem, RPA has been sh own to capture the essential 
physics for doped graphene! 14 * 15 * 21 * contrary to undoped 
graphene, the density of states is finite at the Fermi level 
once doping moves ep away from the Dirac points. How- 
ever, due to the vanishing density of states at the Dirac 
point p 2 ^ the RPA description is incomplete for undoped 
graphene, 23 and it is therefore necessary to consider an- 
other class of diagrams in perturbation theory (such as 
ladder-type vertex corrections) which lead to a plasmon 
resonance below the threshold uo < v^q. For doped 
graphene, the poles of n jRPA (q, uo) define the dispersion 
of a collective plasmon mode, which has the same long- 
wavelength sjq behavior as the plasmon in the 2D elec- 
tron gas studie d in the preceding section [see Eq. (11) 
and Fig. |2(b)| . An approximate dispersion relation of 
the plasmon in a single layer of doped graphene was cal- 
culated in Ref. [12] in the framework of a study of in- 
tercalated graphite and was found to coincide with Eq. 
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FIG. 3: Polarization function of a 2DEG [plots (a) and (b)] 
and graphene [(c) and (d)] for different wave- vectors go- Solid 
lines set for Re II (qo, uo) and dashed lines for Imn°(go,^). 
qo = 0.5/cf in (a) and (c), and qo = 2.5&f in (b) and (d). All 
the plots are done for 5 — > + . 



[TTJ in which the Fermi energy is now related to the uni- 
form density of electrons n e i as ep = vpkp = v-p\/7rn e i. 
Within the previous approximation, the plasmon mode 
enters the inter-band region of the PHES at a momentum 



4/c F 2 + r s 



4r s 



(18) 



where r s = e 2 /e^v-p. Note that r s is density-dependent 
in the 2DEG whereas it is scale- invariant in graphene. 
One important difference with respect to the 2DEG case 
is that only for uo > v^q it is possible to have a solu- 
tion for the plasmon dispersion 1 = v(q)Ke n°(<2, uo), be- 
cause Ren°(q,o;) < for uo < v F q [Fig. |fc)-(d)]. As 
a consequence, the collective plasmon mode in graphene 
at zero magnetic field can only be damped when decay- 
ing into inter-band (and never into intra-band) particle- 
hole excitations. Notice also the difference with respect 
to the 2DEG case, where the plasmon mode is damped 
once its dispersion touches, at some critical wave- vector, 
the boundary of the PHES. The RPA dispersion relation 
of the 2DEG plasmon itself, however, never enters the 
intra-band continuum. In the case of graphene, even if 
the plasmon enters the inter-band region of the PHES at 
a well defined wave-vector approximately given by Eq. 
(18), the mode continues to exist in a rather large re- 
gion of the inter-band continuum! 2 -^ Recently it has been 
argued 25 that, due to the lack of Galilean invariance in 
graphene, exchange interactions, which are not included 
in the RPA, renormalize the plasmon dispersion of doped 
graphene in the long- wavelength limit. This renormaliza- 
tion is due primarily to non-local inter-band exchange in- 
teractions, which reduce the plasmon frequency relative 
to the RPA value. 

We study now the PHES of doped graphene in the 
presence of a strong magnetic field perpendicular to the 



sample. As in the zero magnetic field case, the strong 
contribution to the polarization comes from the diver- 
gence of n°(q, uo) at uo = v^q, as it can be seen in Fig. 



2(c) The effect of the LL wave- function overlap T™ n , (q) 
[Eq. (A7)] is appreciable, in the intra-band region of the 
spectrum in Fig. [4^b), where we see that the intensity 
of the modes is larger near the threshold v^q and prac- 
tically unappreciable near the second boundary of the 
intra-band PHES uo = max(0,^F<7 — 2ep). But in ad- 
dition, Imn°(q, uo) is finite not only in the intra-band 
region, but also in the zones of Fig. |2(c)| with a finite 
weight, as the yellow stripes above and below uo = v^q. 
This form of the PHES is due to both, the LL structure 
of the spectrum and the presence of inter-band exci- 
tations that lead to the mentioned stripes in region II of 
the PHES. The most salient fea ture t hat we find compar- 
ing t he PH ES of a 2DEG [Fig. |l(c)| to that of graphene 
[Fig. 2(c)| in a magnetic field, is that in the former the 
spectrum is composed of horizontal and equidistant non- 
dispersive lines, while in the latter these modes are not 
visible, and the important modes are the diagonal lines 
parallel to the threshold uo = v^q. 

This particular feature of the PHES of graphene in a 
strong magnetic field may be understood in the follow- 
ing manner. Notice first that, in contrast to the 2DEG 
with its equally spaced LLs, the spacing of the relativis- 
tic LLs in graphene decreases at higher energies. In 
a fixed energy window at high energies, there are there- 
fore more possible inter-LL excitations from the level n 
in the band A to n' in the conduction band, of energy 
^nn'-n ~ V%(vf/Ib) ( — A\/ri), than at lower ener- 
gies. Notice further that above an energy of ep also inter- 
band transitions with A = — contribute. Even for small 
values of 5, i.e. in clean samples, neighboring LL tran- 
sitions overlap in energy such that the horizontal lines, 
which dominate the PHES of the 2DEG in a strong mag- 
netic field [Fig. [1(c)] , are blurred. At fixed energy, the 
spectral weight is again not homogeneously distributed. 
This is a consequence of the wave-function overlap be- 
tween the electron and the hole, as we discuss in more 
detail in the following section. 

The above-mentioned dispersive modes in graphene, 
acquire coherence once electron-electron interactions are 
taken into account. This can be seen in Fig. [2(d)] which 
shows the RPA polarization function, and where the dis- 
persive diagonal lines are now clearly distinguishable. 
We will refer to them as linear magnet o-plasmons. In 
the inter-band region of the PHES, the number of lin- 
ear magneto-plasmons depends on the high energy cutoff 
N c . We should keep in mind that the RPA is a good 
approximation for describing the long-wavelength part 
of the spectrum, but fails in reproducing many of the 
physical properties of a system in the short-wavelength 
regime. The dispersion of the collective modes at short- 
wavelength is renormalized when diagrammatic contri- 
butions beyond the RPA are taken into account. In par- 
ticular, the inclusion of ladder diagrams, which account 
for the direct interaction between the electron and hole, 
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FIG. 4: Density plot of Imn°(q, uj) in the low energy sector. 
Plot (a) corresponds to a 2DEG and plot (b) to graphene. 
We have used Nf = 3 and 5 = 0.05 in units of uu c and vf/Ib, 
respectively. Islands (see text) are clearly visible in the oo = uj c 
horizontal line in (a) and in every horizontal line in (b). 



as well as the exchange terms, lead to the excitoni c and 
exchange shifts in the magneto-exciton dispersion ! 3 * 5 * 11 * 



V. STRUCTURE OF THE PARTICLE-HOLE 
EXCITATION SPECTRUM: A WAVE-FUNCTION 
ANALYSIS 

The boundaries of the PHES in a magnetic field may 
be understood by considering an electron-hole pair and 
treating the cyclotron motion of both the electron and the 
hole in a semiclassical limit. The boundaries are related 
to the region in real space where the electron and hole cy- 
clotron orbits may overlap. If we decompose the position 
of an electron r into its cyclotron r\ and guiding center R 
coordinates, r = R + r], the finite overlap of the electron 
and hole orbits implies that n' c — n c < AR < n' c + rj Cl 
where AR = R' — R, AR = |AR| and rj' c and r] c are 
the cyclotron radius of the electron and the hole, respec- 
tively. The latter are given by n c = (\tj\) = ls\/2n + 1 
and r]' c = Ib^/^ti' + 1 in terms of the LL indices n and n! 
of respectively the hole and the electron. As the distance 
between guiding centers AR is related to the electron- 
hole pair momentum q by AR = ql B (see Appendix [b| , 
the momentum is constrained to 



V2n f + 1 - V2n + T < ql B < ^2n' + 1 + V^n + 1. (19) 

As an illustration, the boundaries of the PHES at m = 
n' — n = 1 and for Np — 3 obtained from Eq. jT9] ), 
0.35 < qls < 5.65, coincide with that shown in Fig. Hpi) 
and (b). 

The presence of a set of islands - which is the name 
we give to regions of high spectral weight within each 
electron-hole contribution (n + m, n) to a given horizon- 
tal line m - can be understood from the form of the wave 
functions of the electron and the hole forming the pair. 
In the symmetric gauge, the modulus \^ n ^(x,y)\ of the 
LL wave-function is rotation-invariant, its shape being 
that of n + 1 concentric and equidistant rings (of aver- 
age radius r$ ~ V2£Ib)^ Therefore, one expects that 



m = 1 


























m = 3 





























FIG. 5: Imn°(q,cj m ) for a 2DEG [plots (a), (c) and (e)] and 
for graphene [plots (b), (d) and (f)]. Here uj^ )EG = muj c 
in a 2DEG and u% aph = vfI' 1 ( ^2(N F + m) - V^Nf) in 
graphene. The inset of each figure represents the different 
particle-hole excitations contributing at that energy. As we 
use Nf = 3 for all the plots, there are Nf + 1 = 4 maxima 
in ImII (q, uj m ) corresponding to islands (see text). Notice 
that for graphene, we show the polarizability at th e ener gy of 
the first, second and fourth horizontal line of Fig. 4(b) The 
third line at this filling corresponds to the (Nf + 1, Nf — 1) 
electron-hole pair and presents therefore Nf peaks. 



electron-hole excitations of momentum q will be possible 
whenever there is a finite overlap of the particle and hole 
wave-functions, the guiding centers of which are sepa- 
rated by a distance AR = l B q. If n' is the LL index 
of the electron and n that of the hole, there are n + 1 
substantial overlaps of the rings of the electron and the 
hole wave functions when r)' c — r) c < AR < iq f c + r\ c . This 
will lead to a division of the contribution (n 7 ,n) to the 
m-th horizontal line of the PHES into n + 1 regions of 
preferred momenta or islands. 

In order to understand the effect of these overlap func- 
tions in more detail, we first consider the line m = 1 of 
lowest energy. Both in the 2DEG and in graphene, the 
only possible inter-LL excitation that contributes to the 
formation of this energy line involves a hole in the LL 
Nf and an electron in Np + 1. T his is schematically rep- 
resented in the inset of Fig. |5(a)f (b), where we represent 
the unique electron-hole transition contributing to the 
first horizontal line in each case. As argued above, there 
are Nf + 1 zones of preferred momenta due to the over- 
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lap between the wave functions of the electron and the 
hole. In Fig. [HJa) and (b), we have plotted ImII (q, uj m ) 
for the 2DEG and graphene, respectively, at the energy 
corresponding to the first horizontal line, with Np = 3. 
One obtains indeed four peaks which yield the islands 
observed in the low-energy zoom of the PHES (see Fig. 
[4] where we have also chosen Np = 3). Notice, however, 
that the last island, though present, is hardly visible in 
graphene due to the suppression of backscattering. 



At larger values of m, there is an essential difference 
between the 2DEG and graphene. Due to the equidis- 
tant level spacing in LL quantization for non-relativistic 
electrons in the 2DEG, the horizontal lines in the PHES 
occur at uo m = mco c . If m > 1, there are min(m, Np) dif- 
ferent inter-LL transitions (n + m, n) that contribute to 
the spectral weight of the m-th horizontal line [see the in- 
set of Fig. 5(c)| and (e) which indicates the electron-hole 
transitions contributing to the second and third horizon- 
tal line in a 2DEG]. But all these transitions have differ- 
ent preferred momenta because of different electron-hole 
overlaps. The spectral weight is therefore a superposi- 
tion of these different overlap functions [see Fig. [5|c), (e), 
where we plot ImII (q, uo m ) for m = 2, 3] and the islands 
are no longer well defined in the horizontal direction, as 
it may be seen in the second and third horizontal line of 
Fig. |4(a)| which have lost the dashed structure of the first 
line. The situation is remarkably different in graphene, 
where the LL spacing is not constant and where a par- 
ticular horizontal line is due to the inter-LL transition 
with energy uj^ m = \/2(^f /lB)W n + 171 ~ ^V™) an d 
therefore not only determined by the LL-index separa- 
tion m. Apart from very rare events in the high-energy 
regime where two inter-LL transitions (ni + mi, n\) and 
(n 2 + m 2 ,n 2 ) may coincide in energy cj* 1>mi = ^ )TO2 , 
each horizontal line therefore consists of a single inter- 
LL transition and has Np + 1 well-separated peaks in 
ImII (q, cc^ m ), as we have shown in Fig. |H|d), (f) for 
(N F + m,m), with m = 2,3. In contrast to the 2DEG, 
the islands remain thus well separated in the horizontal 
direction, whereas they overlap strongly in the vertical 
direction (i.e. in energy) due to the decreasing level spac- 
ing at higher energies and the large number of inter-LL 
transitions in a fixed energy window [Fig. |4(b)| , as we 
have discussed in the last section. As a consequence, the 
most prominent modes in graphene in a strong magnetic 
field are diagonal lines, parallel to uo = v^q, whereas those 
in the 2DEG remain horizontal. Electron-electron inter- 
actions turn these lines of large spectral weight into co- 
herent modes: magneto-excitons in the 2DEG and linear 
magneto-plasmons in the case of graphene, in addition to 
the upper hybrid mode that reveals itself in the formerly 
forbidden energy region of the PHES of non-interacting 
particles. 



VI. STATIC SCREENING 

In this section we study the properties of n°(q) = 
n°(q, uo = 0) in the static limit, for which the polarization 
is entirely real. The polarizab ility of graphene is shown 
in Fig. ga) and (b) for B = d 12 ' 14 ' 15 ' and B ^ 0, respec- 
tively. For comparison, we also show the corresponding 
polarizability of a standard 2DEG 10 in the absence [Figs. 
|6(c)| and in the presence [Fig. |6(d)| of a magnetic field. 
In order to compare the B 7^ to the B = polariz- 
ability, we have chosen a Fermi wave-vector kp for the 
B = case equal to \/2Np + 1/Ib, which corresponds to 
the same carrier density as a graphene layer in a mag- 
netic field with all LLs filled up to the A^-th level of the 
conduction band. 

One first notices a difference in the B = static po- 
larizabilities between graphene and the 2DEG. Although 
the static polarizability remains constant and equal to 
the electronic density of states p(ep)^m both cases up 
to a wave vector 2k p, there are two contributions for 
graphene that stem from intra-band and inter-band ex- 
citations, respectively. Whereas the polarizability due 
to intra-band excitations in graphene [red dotted line 
in Fig. [6^a)] decreases linearly in q, due to the elec- 
trons' chirality (17) and the absence of backscattering, 



the intra-band contributions yield a linearly increasing 
polarizability. Beyond q = 2kp, there are no possible 
zero-energy particle-hole excitations in the intra-band re- 
gion, and the associated polarizability there fore t ends to 
zero. This is also the case in the 2DEG [Fig. 6(c) , where 
there are only intra-band excitations. In graphene, how- 
ever, inter-band excitations still yield a linearly increas- 
ing contribution to the total polarizability, which then 
asymptotically approaches the inter-band polarizability 
[blue dashed line in Fig. [6^a)]. 

Qualitatively, one finds a similar behavior for the static 
polarizability at 5 / except in the small-g limit. 
Whereas the static polarizability at B = remains con- 
stant and coincides with the density of states at the Fermi 
energy, it tends to zero as n°(q — > 0) ex q 2 for B 7^ 0P^ 
This is due to the fact that the main contribution to the 
polarizability comes from q = excitations in the vicin- 
ity of the Fermi energy ep. Contrary to the B = case, 
where there are q = excitations the energy of which 
tends to zero, ep lies now in the cyclotron gap between 
the highest occupied LL Np and the lowest unoccupied 
one Np + 1. This energy gap must be overcome by q = 
excitations, such that its spectral weight tends to zero 
then. Indeed, the static polarizability also coincides with 
the density of states at the Fermi energy because the 
latter vanishes for B ^ when ep lies in the gap. 

Furthermore, one notices the oscillatory behavior of 
the static polarizability, both for graphene and the 
2DEG, below 2kp. These oscillations are again due to 
the wave-function overlap between the electron and the 
hole, and one obtains Np + 1 maxima. Since the main 
contribution to the polarizability at small wave vectors 
comes from excitations in the vicinity of ep, the oscil- 
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FIG. 6: (Color online) (a) Vacuum (inter-band) contribution 
(blue dashed line), intra-band contribution (red dotted line) 
and total static polarization function (black thick line) of 
graphene for B = (see ReffTB]). (b) Same as (a) but for 
B^0. We have used Nf — 3 and N c = 350. For consistency, 
the Fermi momentum used in (a) is kF = V^Nf + 1/Ib = 
V^/Ib- (c)-(d) Same as (a)-(b) respectively, but for a stan- 
dard 2DEG. 



lations are dominated by the (Np + l,iVp) intra-band 
transition in graphene, as one may also see from the red 
dotted line in Fig. [6^b), which represents the intra-band 
contribution to the polarizability. At large values of the 
wave vector q the static polarizability is, as in the B = 
case, dominated by inter-band excitations the discrete 
nature of which is less important than in the small-q 
limit. The linear increase therefore coincides with the 
B = result. 

The static polarizability is a useful quantity for the cal- 
culation of the screening properties of electrons. Screen- 
ing, e.g. of the Coulomb interaction potential between 
the electrons or the potential of a charged impurity, 
is indeed determined by the (static) dielectric function 
e RPA (q,u = 0) = 1 - v(q)U°(q). At zero field, the long 
wavelength limit is similar in the two cases: e RPA (q) « 
1 + Qtf/q^ where qpF = 2i\e 2 p(ep) j e\> is the Thomas- 
Fermi wave-vector. Note however that qpp is density 
independent in the 2DEG whereas it scales as kp in 
graphene. Therefore, the dielectric function diverges as 
e ~ q TF /q — > oo in the two cases when q — > 0. However, 
the kp dependence in the numerator of e RPA (q — > 0) 
in graphene points out the absence of screening in un- 
doped graphene (kp = 0) at long distances. In the short 
wavelength region q ^> 2kp, e RPA (q) tends to 1 in a 
2DEG, whereas for graphene, e RPA (q) — > l+7rr s /2. This 
extra contribution of 7rr s /2 to the dielectric function of 
graphene at large wave vectors is due to the linear growth 
of the polarizability at q > 2kp and is therefore related 
to virtual inter-band particle- hole excitations.^ In sum- 
mary, at short wavelengths, a 2DEG does not screen at all 




FIG. 7: Static (oo = 0) dielectric function for graphene (a) 
and a 2DEG (b) in a magnetic field computed in the RPA. 
We have used Nf — 1,2,3, corresponding to the blue, red 
and green lines, respectively. 



[e —> 1), whereas (doped or undoped) graphene screens 
as a dielectric (e —> 1 + 7rr s /2 > 1) thanks to its filled 
valence band. 

The situation is different in the presence of a magnetic 
field. In Fig. [7ja) and (b) we have plotted the static 
dielectric function for graphene and the 2DEG, respec- 
tively, in a magnetic field (see also Ref. 4 26]). Notice that 
at long wavelengths, e RPA (q) — 1 ex q in the 2DEGf^ as 
well as in graphene. In fact, in the limit of Np ^> 1, 



„RPA 



(?)-!« r s Ny 2 ql L 



(20) 



as q — » 0, this limit being valid for both, a 2DEG and 
graphene. The difference stems again in the density de- 



pendence of r s in the two cases: because r s 



-1/2 



in a 2DEG, e RPA grows linearly with Np in this case. 
However, r s is density-independent in graphene, lead- 

3/2 

ing to a dielectric function proportional to . Fur- 
thermore the maximum of the dielectric function be- 
haves as £ max e RPA (q ~ l/kpl 2 B ) ~ r s N F . There- 
fore £ max cx \/Np in a 2DEG whereas £ max oc Np in 
graphene. This different behavior is reflected in Fig. 7(a)| 
and (b). As we see, there is a considerable increase of the 
static dielectric function of graphene as we increase iVp, 
as compared to the 2DEG. This is due to the relativis- 
ts LL quantization of graphene, which leads to an in- 
creasing of the quantum effects (virtual inter-level tran- 
sitions) as the separation between levels becomes nar- 
rower. On the other hand, in both graphene and the 
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2DEG, e RPA (q,u = 0) 1 as q 0, which implies 
that there is no screening at long distances (as in the 
undoped zero field case, as discussed above). The short 
wavelength behavior of the dielectric function in a mag- 
netic field is, however, similar in both the 2DEG and 
graphene to their respective zero field limits. Therefore, 
the short wavelength decay of the effective interaction 
in graphene, due to inter-band polarization effects, leads 
to a screening similar to that of an insulator, while the 
intra-band processes leads to a metallic-like screening. 



VII. CONCLUSIONS AND OUTLOOK 

In conclusion, we have compared in detail the polariz- 
ability of doped graphene with and without a strong mag- 
netic field to that of the 2DEG. In the absence of a mag- 
netic field, the main difference arises from the presence 
of two different regions of non-vanishing spectral weight 
in the PHES ImII(q, uj) of graphene. These two regions 
represent contributions from intra- and inter-band exci- 
tations, respectively, whereas in the 2DEG with only one 
parabolic band, there is only one region. Furthermore, 
the chirality of electrons in graphene suppresses backscat- 
tering such that the spectral weight is centered around 
the main diagonal of the PHES at uo — v^q, whereas it is 
more or less homogeneously distributed over the particle- 
hole continuum in the 2DEG. Electron-electron interac- 
tions yield in both cases a plasmon mode that disperses 
as J~q. 

In the presence of a strong magnetic field, the plasmon 
mode evolves into the upper hybrid mode which is gapped 
at zero energy. In the 2DEG, this gap is given by the cy- 
clotron frequency uo c1 whereas in graphene it is eBv^/ep 
and thus depends on the Fermi energy ep- The most 
salient difference between graphene and the 2DEG are 
the additional modes that occur in the interacting system 
in the parts of the PHES that correspond to the zero- 
field particle- hole continuum. In the 2DEG, the spec- 
tral weight is concentrated along equidistant horizontal 
lines, due to the equidistant LL spacing, and the result- 
ing magnet o-excitons are therefore weakly dispersing. In 
contrast to these rather well studied modes in the 2DEG, 
one finds linear magneto-plasmons in graphene that dis- 
perse roughly parallel to the central diagonal uo — v^q in 
the PHES. These modes are a consequence of a different 
organization of the regions of highest spectral weight in 
graphene as compared to the 2DEG. The energy levels in 
graphene are no longer equally spaced as a consequence of 
relativistic LL quantization, and the energies of inter-LL 
transitions are more densely packed than in the 2DEG, 
especially at higher energies. Even a small level broad- 
ening due to impurities therefore leads to an overlap in 
energy of the inter-LL transitions. Furthermore, spectral 
weight is highly modulated at a fixed energy due to the 
wave-function overlap between the electron and the hole 
involved in the excitation. 

We have finally discussed the static polarizability and 



the dielectric function that describe the screening prop- 
erties of the system. In contrast to the 2DEG, where 
the polarizability tends to zero at wave vectors larger 
than 2/cf, it increases linearly in graphene, due to the 
increasing relevance of inter-band excitations. Whereas 
the small-*? behavior of the dielectric constant is similar 
in graphene and the 2DEG, a calculation within the RPA 
indicates that it tends to a constant different from one in 
the large-g limit for graphene. 

As for an experimental confirmation of the partic- 
ular high-field collective excitations in graphene dis- 
cussed above, one may first think of magneto-optical 
experiments. Transmission spectroscopy has indeed re- 
vealed the characteristic V Bn behavior of the graphene 
LLs in epitaxial and exfoliated graphenef^ in agree- 
ment with theoretical expecations ! 30 * 31 * 3 -^ Similarly, Ra- 
man spectroscopy has been successfully applied to 
graphene in the absence^^ and in the presence of a 
magnetic field l 3 -^ Thus the particular electron-phonon 
interactio n ^ 36 * 37 * a nd the theoretically studied magneto- 
phonon resonanc e 138 * 39 * 40 * could be confirmed. However, 
these techniques are restricted to zero wave- vector excita- 
tions, whereas electron-electron interactions and the re- 
sulting collective excitations are more prominent at non- 
zero values of the wave vector. In order to probe the 
excitation spectrum at non-zero values of the wave vec- 
tor, inelastic light scattering may be a promissing tech- 
nique that has been successfully used t o stud y collective 
quantum-Hall excitations in the 2DEG. E ^1 
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APPENDIX A: CALCULATION OF THE 
POLARIZATION FUNCTION 

The electronic wave function in graphene in a magnetic 
field can be expressed as a four component spinor, 



*+(r,t) =e iK,] 



tf_(r,t) = e- iK,] 




/ 




(Al) 



(A2) 



with the components 

$ Ca (r,t) ^ ^2( T \^<;a;\n£)c<; ; \n£(t) 



(A3) 



A=± n.t 
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where c^ ; \ n £(t) is the annihilation operator of an electron 
in the state |n, £) of the A band with valley index £. The K 
valley (£ = +) part of the single particle Green's function 
in reciprocal space can be written as 



G°A<h«>) = 



g \,aa(^u) 



G \,ba(<±^) G° +BB (q,uo) 
with matrix elements 



(A4) 



EE 

A n 



olol' ;An 



(k + CK) 



UO ' 



A£ n + Wsgn(A^ n ) ' 



(A5) 



where A£ n = Ae n — ep is the energy difference between 
the LL and the Fermi energy ep — e at f , which we choose 
in the conduction band (A = +). Furthermore, S is a 
positive infinitesimal in the clean limit and the matrix 



f\ n (<l) has been derived in Ref. |8| The expression for 
the Green's function in the K' valley (£ = — ) may be 
obtained with the help of G®. aa , = G°_£. a , a , and the non- 
interacting polarization operator II (q, uo) can be calcu- 
lated from Eq. Q. Taking into account both valleys, 
the particle-hole polarization reads 



in°(q,w) = i [11° (q,w) + 11° (q,w)] = 2ffl° (q.o;), 

(A6) 

where the last step indicates that one obtains equal con- 
tributions from both valleys. We may therefore restrict 
the calculation to only one valley (£ = +), and take into 
account the twofold valley degeneracy by a simple factor 
9v = 2. 

The integration over the frequency integral yields then 
Eq. where the functions T nn , (q) are given by 



-AA' 



-l 2 B q 2 /2 

2ttZ| 



I 2 a 2 
i B q 



Al* 1*/^ 



D! 



(n>-l)! 



r n>- 



n< 
1 



r 



/ 2 tf 2 



A/r)* r)* 
Z n Z n' 




rn>- 



I 2 a 2 



iAZl 



We define^ 



n 



AA' 

+ ,nn' 



where o; + — > 
— — z5 and 



A£n - A'f n / + + W 

indicates the replacement uo + z5 



a;") 
(A8) 



However, due to the fact that the separation between 
LL in graphene decreases with n, it is always possible to 
have semiquantitative good results from smaller values 
of N c . 



n-l 



A' n'=0 

+ E E n + ^(q,c) 



A' n'=n+l 

+ n^- n A „(q,uO 



(A9) 



which verify II* n (q, uo) = — 11^ (q, uo). The vacuum 
polarization, which accounts for the inter-band processes, 
is defined as 



ffnq^) = -E n +7n(q^) 



(A10) 



where N c is a cutoff. Taking into account that, al- 
ready in the absence of magnetic field, the validity of 
the continuum approximation is up to A ~ t, then 
6_/v c = (v F /l B )V27Tc ~ which leads to iV c - 10 4 /B[T], 
which is very high even for strong magnetic fields. 



APPENDIX B: ELECTRON-HOLE PAIR 
MOMENTUM 

In this appendix, we relate the momentum q of an 
electron-hole pair to the distance between the guiding 
centers AR of the electron and the hole. In classical 
mechanics, the cyclotron motion of an electron leads to 
7r = eB x 77. where 7r is the gauge- invariant momentum 
and rj is the cyclotron coordinate. From Newton's equa- 
tion with Lorentz's force, it is obvious that the quantity 
K = 7T — eB x r is a constant of the motion, where r is the 
electron position. This constant of the motion is usually 
called the pseudo-momentum or generator of magnetic 
translations.^ Defining the guiding center coordinate as 
R = r — 77, the pseudo-momentum reads K = eR x B, 
which actually shows that, apart from a conversion fac- 
tor eB = 1/1%, the guiding center coordinate and the 
pseudo-momentum correspond to the same constant of 
the motion. 

Now consider an electron- hole pair, where the elec- 
tron has a pseudo-momentum K' and the hole a pseudo- 
momentum — K (corresponding to a removed electron of 
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pseudo- momentum K). The pair has a momentum 

q = K' - K = eAR x B (Bl) 



and therefore AR = ql 2 B , which is the sought after rela- 
tion. 
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